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We represent by an electrical model the imagine} of a one-dimensional 
coherently illuminated and diffusely reflecting surface by an optical system 
with a rectangular aperture. We then obtain the statistical properties of 
the image intensity from the statistical properties of the square of the en- 
velope of a narrowband noise signal in the electrical model. The analysis 
is simple because use can be made of results known in communication 
theory. The results agree with those obtained in a direct way. 

I. INTRODUCTION 

The speckle pattern in the image of a coherently illuminated and 
diffusely reflecting object has been analyzed by Enloe 1 . Enloe's results 
show a remarkable similarity to some results occurring in the theory of 
narrowband noise (See Ref. 2, pp. 397 ff.). This similarity prompts the 
question whether Enloe's results can be derived by the use of an elec- 
trical model analog}' involving narrowband noise.* In this paper we will 
show that this is indeed possible for the special case of a one-dimensional 
subject and an optical system with a rectangular aperture. The analysis 
is simple because we can use results known in communication theory 
and the results agree with those of Enloe 1 . 

II. THE OPTICAL MODEL 

The optical model is shown in Fig. 1. In plane Pi there is a coherently 
illuminated row of scattercrs which scatter with random phase (a one- 
dimensional diffusely reflecting surface). At a distance d in plane P 2 
there is a lens. In front of the lens there is a rectangular aperture with 
an amplitude transmission H(x 2 , */ 2 ). The distance d is large compared to 
the focal length / of the lens and therefore to a good degree of approxima- 



* Such an analogy was already suggested by Rigden and Gordon. 3 
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Fig. 1 — Optical model. 

tion the image of the scatterers in plane Pi is situated in plane P 3 at a 
distance / from the lens. Enloe computed the intensity as well as the 
autocorrelation of the intensity in the image plane. 1 

As a help in understanding the following electrical model we make 
the following comments regarding the optical model. Suppose that the 
instantaneous value of the electric field e(x, , t) in the object plane Pi 
is given by 



e(aj, , = o^l exp (-2tt> o 0- 
e(xi , t) is zero for y x ^ 0. The time-independent phasor is 

a \u) = a(a>) ' 

where we have introduced the spatial frequency coordinate 

Xt 



(1) 



(2) 



(3) 



(We write a( ) as a function of the spatial frequency xj\d and not of 
Xi in order to simplify the following computation.) For later use we also 
introduce 



Hi. 



(4) 



Since the lens is situated in the far field of the object, the phasor A (x 2 , yi) 
in plane P 2 in front of the aperture is the Fourier transform of a(cti), 



A(x 2 ,y 2 )=j I a(«!) exp (2Trja 1 x 2 ) exp (2tj/3,t/ 2 ) da x d/3, . 



(5) 
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Since o(«i) is zero for 0, ^ the above two-dimensional Fourier trans- 
form operation reduces to the one-dimensional Fourier transform 
operation 



A(x 2 , y 2 ) = / a(«i) exp (2TJa l x 2 ) da, 

J —CD 



(6) 



It is seen that A(x 2 , y 2 ) is a function of x 2 only and we will therefore 
write A(x 2 ). The phasor behind the aperture, B(x 2 , y 2 ) is obtained as 
the product of A(x 2 ) and the aperture transmission function H(x 2 , y 2 ) 

B(x 2 , y 2 ) = A{x 2 ) -H(x 2 , y 2 ). (7) 

As was mentioned we assume a rectangular aperture function. We 
assume that H(x 2 , y 2 ) can be written as 

H(x 2 ,y 2 ) = H^xJ-HM. (8) 

Since our object is one-dimensional we are only interested in the image 
intensity at y 3 = 0. The image intensity at y 3 = is not influenced by 
H 2 (y 2 ) except for a factor which remains constant over all x 3 . For 
equation (7) we can therefore write 



B(x 2 ) = AixJ-H^x,) 



(0) 



with the understanding that B(x 2 ) does vary in the y 2 direction but that 
this variation is of no interest to us. So much for the optical model. 
We will now present the electrical model. 



III. THE ELECTRICAL MODEL 

The electrical model is shown in Fig. 2. The electric field in the object 
plane is scanned by a detector whose output voltage is proportional 
to the instantaneous value of the electric field in the object plane. (In 
the optical region there are no such detectors available. This should 
not present any conceptual difficulties since we can scale up our model to 
a longer wavelength where there are such detectors.) 




Fig. 2 — Electrical narrowband noise model. 
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The electric field e(xi , t) in plane Pi is given by equation (1) and the 
output voltage v,(/) of the scanning detector is 



v ( (t) = Cial^j exp (—2irjv t). 



Vi is the scanning velocity and we have 

x, = v,t. 



(10) 



(ID 



c x denotes a constant of proportionality which is of no interest. The 
Fourier transform of i/,-(<) of equation (10), FT[Vi(t)] is given by 



FT[v { (l)] = Cl — A 



£(? -v )J 



= c,A 



Mfr-4 



(12) 



.-l(j') is the Fourier transform of a(t). c 2 is a constant of proportionality. 
The Fourier transform or spectrum of v t (t) is centered at v as shown in 
Fig. 3. In the optical model the spectrum A(x 2 ) is multiplied by the 
aperture transmission function H l (x 2 ) in plane P 2 . See equation (9). 
In order to simulate this in our electrical model we have to pass the 
voltage Vi(t) through a temporal frequency filter with the frequency 
response Hj[()d/Vi)(v — v )]. The spectrum B[(\d/Vi)(v — v )] at the 
output of the filter is then given by 



6r<— >] 



B — 



= A 



["fr -"•>]*■!)?<' -4 < i3 > 



The filter i7i[(\d/vi)(y — v )] is shown in the electrical model of Fig. 2 
and its frequency response is sketched in Fig. 3. (The frequency response 
ffi[(\d/vi)(v — v )] in the electrical model corresponding to the rectangu- 
lar aperture transmission function of the optical system can only be 
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Fig. 3 — Spectrum of the electrical signal and frequency response of the filter. 
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realized with a time lag, but this need not disturb us.) At the output of 
the filter we have the voltage v n (t) 



,M^1 e>:!> -- 



(14) 



where 6(0 is the Fourier transform of B(v). c 3 again is a constant of 
proportionality of no interest. The output device scans the image plane 
P a with a velocity v 2 and 



(15) 



because in the optical model the image is demagnified by a factor cl/f 
with respect to the object. The inversion of the optical image with re- 
spect to the object is accounted for by inverting the coordinate system 
in the image plane P 3 , see Fig. 1. We can now write for equation (14) 



r„(D = *&(*J 



exp (— 2n>„0 



= C 8 b(^) exp(-2irju t). 



(16) 



In the optical model we detect the intensity in the image plane, that is, 
the square of the absolute value of the phasor. Therefore in our electrical 
model the output device displays | 6(.r 3 /X/) | . 

The voltage v,(t) at the output of the filter is narrowband compared 
to j*„ . We now consider !•„(/) as a narrowband noise voltage. Another 
equivalent representation is (see, for example, Ref. 2, p. 397) 

,.„(/.) = N r (0 cos LW + N.(t) sin 2W, (17) 



where 



N f (t) = Re 
N.(t) = Im 



My 

Mi 



(18) 



Re [ ] means "real part of" and Im [ ] means "imaginary part of". 
We now assume that the filter is so narrowband that its impulse response 
is much wider than the time over which y,(i) shows any appreciable 
correlation. (This assumption is the equivalent of Enloe's assumption 
that the optical spread function is much wider than the average distance 
between the images of scatterers with independent phase.) If this 
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assumption holds, the voltage v (t) results from many independent 
values of v { (t). According to the central limit theorem v (t) and therefore 
N e (t) and N,(t) show a Gaussian distribution. The voltage v (t) therefore 
has the statistical properties of narrowband Gaussian noise which are 
discussed for example in Ref . 2, pp. 397 ff. As was mentioned, the square 
of the envelope corresponds to the optical intensity. The envelope E(t) 
is obtained as 

E(t) = [N 2 c (t) + Nl(fff 

,.\i (19) 



({*• [»&)]}' + {- M 



= lt Re 

According to Ref. 2, equation (9.18), E(t) has a Rayleigh probability 
distribution density W(E) 

_ E exp (-g/2,) 

The average value of the square of the envelope (E 2 ) av which is equal to 
the average value of the optical intensity (I) av is given by [Ref. 2, Eq. 
(9.21a)] 

{E\. = (I) av = 2+. (21) 

We are now also interested in the autocorrelation R(s) of the intensity / 

R(s) = (I(x 3 )I(x 3 + s))„ (22) 

= (E\x 3 )E\x z + 8))„ . 

According to Ref. 2, equation (9.24), R(s) is given by the following 
expression 

R(s) = 4* 2 [1 + k 2 (s)] (23) 

= </ 2 )av[l + «(«)], 

where we have used equation (21). As explained in Ref. 2, [equations 
(9.10b and 9.12b)] k (s) is the autocorrelation of the spread function 
corresponding to a filter with the frequency response H^Xd/v^p). 
This is the frequency response of the filter in our electrical model (see 
Fig. 2), but centered at zero frequency. The frequency response fl r 1 ((Xd/»i)i») 
is shown in Fig. 4. Jc (s) is normalized to one at s = 0. (The above dis- 
cussion holds when the frequency response of the filter is symmetrical 
about the origin as is true for our case.) The spread function is the 
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Fig. 4 — Frequency response of the filter HA(\d/vi)v]. 
Fourier transform of the frequency response and so we have 



«M£' 



-«*Q 



(24) 



where H^v) and h(t) are Fourier transform pairs. c 4 is a constant of no 
interest. Therefore we obtain for k (s) 



™-3sC*®fii+% 



p(0) 

fV,T 



di 



(25) 



p(0) p(0) ' 



where 



p(tt) = | /i *(/)/«(/ + u) r//. 
Using equations (23) and (25) we finally obtain for R(s) 



(26) 



H(8) = (I 2 ), 



1 + 



P 2 (0) J 



(27) 



This last equation is the same as Enloe's equation (14) if the last term 
in that equation is disregarded. The last term in Enloe's equation (14) 
vanishes if the spread function of the filter is much broader than the 
average distance between images of scatterers with independent phase. 
The power spectrum follows from the autocorrelation function in the 
same way as in Enloe's analysis and this derivation will not be repeated 
here. 

In summary: We have derived the statistical properties of the intensity 
in the image of a one-dimensional coherently illuminated and diffusely 
reflecting subject with the help of a narrowband noise model. Use was 
made of results known in communications theory. 
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